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ABSTRACT 


In  an  effort  to  solve  the  forward  propagation  problem  associated  with  the  1992  Barents  Sea 
Polar  Front  Experiment  Tomography  Test,  the  transmission  of  a  224  Hz  pulse  signal  from  a  near 
bottom  sound  source  to  a  vertical  hydrophone  array  was  simulated  based  on  three-dimensional  ray 
theory.  Through  numerical  raytracing,  followed  by  eigenray  searches  and  estimations  of  ray 
amplitudes,  phases  and  travel  times,  the  arrival  structure  as  a  function  of  time  and  elevation  angle 
was  constructed.  The  simulation  was  performed  for  both  a  two-dimensional  and  three-dimensional 
modeled  ocean  enviromnent  in  order  to  examine  the  significance  of  three-dimensional  effects.  The 
predicted  arrival  structures  compare  well  with  the  observed  data.  Three-dimensional  effects  proved 
to  be  significant  only  for  the  latest  arrivals. 
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L  INTRODUCTION 


A.  COASTAL  ACOUSTIC  TOMOGRAPHY 

Ocean  tomography  is  an  acoustic  method  fen*  monitoring  the  ocean  interior.  It  uses  low 
frequency  sound  energy  to  collect  information  on  changes  in  the  ocean  temperature  and 
current  stractuie  while  propagating  through  a  water  column.  The  data  collected  are  in  the 
form  of  changes  in  pulse  travel  times.  Inverse  techniques  are  applied  to  this  data  to  create 
a  best  estimate  of  the  unknown  ocean  structure  (Munk  and  Wunsch,  1979). 

A  tomogrs^hy  problem  consists  of  two  parts:  the  forward  problem,  which  establishes 
the  physical  relationship  between  the  data  and  the  unknown  structure;  and  the  inverse 
problem,  which  deals  widi  the  reconstruction  of  the  unknown  environment  based  on  the 
physical  relationships  established  in  the  forward  problem.  The  feasibility  of  using  acoustic 
tomography  to  monitor  the  ocean  interior  depends  upon  the  following  requirements  (Munk 
and  Wunsch,  1979): 

1.  Stability 

2.  Resolvability 

3.  Identiflability 

4.  Adequate  Signal-to-Noise  Ratio 

Stable  and  resolvable  multipath  arrivals  are  necessary  over  successive  transmissions  in 
order  to  derive  a  useful  time  series  of  acoustic  travel  times.  Additionally,  identiflability, 
i.e.,  association  of  modeled  data  to  observed  data,  is  needed  in  order  to  detennine  the 
geometry  of  the  acoustic  paths.  Signal  to  noise  ratio  determines  the  accuracy  of  the  travel 
time  measurements  which,  in  turn,  determines  the  reliability  of  the  acoustic  images  of  the 
ocean  interior. 

In  the  past,  ocean  tomography  was  mainly  used  for  deep  ocean  monitoring.  The 
ad{q)tation  of  the  deep  water  technique  to  shallow  coastal  waters  is  the  main  thrust  of  a 
recent  collaboration  between  scientists  and  engineers  at  the  Naval  Postgraduate  School  and 
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Woods  Hole  Oceanographic  Institution  (Chiu,  et  al.,  1993(a)  and  Miller,  et  al.,  1993). 
Relative  to  deep  water  propagation,  the  multipath  arrivals  from  a  pulse  excitation  in 
shallow  water  tend  to  overlap  in  time,  making  it  difficult  to  resolve  them  individually  with 
the  typical  receiver  configurations  used  for  open  ocean  applications.  In  the  1992  Barents 
Sea  Coastal  Tomography  Test,  Chiu,  et  al.  (1993(a))  and  Miller,  et  aL  (1993)  successfully 
used  a  vertical  line  array  receiver  and  beamfcaming  techniques  to  overcome  this  difficulty. 
Additionally,  Chiu,  et  al.  (1993(a))  developed  new  hybrid  ray-noode  inverse  technique  to 
provide  optimal  sound  speed  maps  with  minimum  mean  square  errors.  As  a  result,  the 
Barents  Sea  Polar  Front  Experiment  has  successfully  addressed  the  requirements  for 
coastal  applications  of  acoustic  tomography. 

B.  BARENTS  SEA  POLAR  FRONT  EXPERIMENT 

The  Barents  Sea  Polar  Front  Experiment  (BSPFEX)  was  conducted  in  August  1992  by 
researchers  from  the  Naval  Postgraduate  School  (NPS),  Woods  Hole  Oceanographic 
Institution  (WHOI),  and  Science  Applications  International  Ccnporation  (SAIC),  (Barents 
Sea  Polar  Front  Group,  1992).  BSPFEX,  conducted  within  an  80  X  90  km  rectangle 
centered  50  km  east  of  Bear  Island,  utilized  both  acoustical  tomography  and  traditional 
oceanographic  techniques  to  study  the  Barents  Sea  Polar  Front 

The  Barents  Sea  is  bordered  to  the  nmth  by  Spitzbergen  and  Franz  Joseph  Land,  to 
the  south  by  Scandinavia  and  Russia,  to  the  east  by  Novaya  Zemlya  and  is  open  to  the  west 
Currents  in  the  region  transport  water  of  both  Arctic  and  Atlantic  origin  that  converge  at 
the  Barents  Sea  Polar  Front  The  colder  Arctic  water  flows  from  the  north,  while  the 
warmer  Atlantic  water  flows  in  from  the  southwest  These  water  masses  have  vastly 
different  tonperature  and  salinity  characteristics  (see  Table  1),  and  the  confluence  of  these 
flows  results  in  complex  vertical  and  horizontal  structure  within  the  frontal  zone.  Hgure  1 
depicts  Loeng’s  (1991)  model  of  the  regions  occupied  by  the  major  water  masses.  The 
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Figure  1.  Barents  Sea  water  masses  (Loeng,  1991) 
abbreviations  are  as  follows:  Arctic  Water  (AW),  North  Atlantic  Water  (NAW),  Svalbard 
Bank  Water(SBW),  Barents  Sea  Water  (BSW),  and  Bottom  Water  (BW).  The  area 
occupied  by  BSPFEX  is  also  indicated  in  the  figure  by  a  rectangle. 

The  physical  oceanogn^hic  characteristics  of  the  experimental  region  have  a 
significant  impact  on  sound  propagation.  The  region  is  shallow,  with  an  average  depth  of 
224  m.  The  bottom,  q)preciably  rough  in  scrnie  areas,  consists  of  thin  sediments  overlying 
hard  rock.  The  range  dependence  of  the  bathymetry  and  sound  speed,  coupled  with  the 
shallowness  of  the  water  depth  can  result  in  complex  sound  refraction  and  boundary 
interacticHis  (EUiot,  1992).  The  effects  of  these  environmental  factors  on  sound  propagation 
will  be  direussed  further  in  Qiapter  n. 
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Table  1:  BARENTS  SEA  WATER  MASSES  (Leong,  1991) 

Mtyor  l^rpes 

T(®C) 

S(psu) 

Arctic  Water  (AW) 

North  Atlantic  Water  (NAW) 

<0.5 

-0.5  to  2.0 

34.2  -  34.0 
34.8  -  35.0 

Local  Variants 

T(®C) 

S(psu) 

Svalbard  Bank  Water  (SBW) 
Bottom  Water  (BW) 

Barents  Sea  Water  ^SW) 
Polar  Front  Water  (PW) 

1 .0  to  3.0 
<-1.5 
-1.5  to  2.0 
-0.5  to  2.0 

<34.4 

>35.0 

34.7  -  35.0 

34.8  -  35.0 

The  Barents  Sea  Polar  Front  Experiment  consisted  of  three  major  components,  a  physical 
oceanographic  investigation  of  the  polar  front,  a  coastal  tomography  test,  and  an  acoustic 
propagation  study.  The  scientific  objectives  of  the  coastal  tomography  test  are  summarized  below 
(Chiu,etal.,  1993): 

1. To  study  the  dynamics  of  the  front  near  Bear  Island  by  mapping  its  spatial  and  temporal 
variability. 

2.  To  assess  the  performance  of  a  new  coastal  acoustic  tomography  system  in  mapping 
frontal  features. 

Figure  2  shows  the  experimental  configuration  for  the  test  The  tomography  system  included 
a  224  Hz  source  (A),  a  vertical  array  of  sixteen  hydrophones  (B),  and  two  400  Hz  transceivers  (C 
and  D).  The  placements  of  the  tomographic  elements  are  summarized  in  Table  2.  This  forward 
modeling  work  focus  solely  upon  the  propagation  path  between  the  224  Hz  soiurce  (A)  and  the 
vertical  hydrophone  array  (B). 

The  224  Hz  source  was  moored  20  m  above  the  sea  floor.  The  vertical  array  consisted  of  16 
evenly  spaced  hydrophones  at  10  m  apart  with  the  bottom  hydrophone  being  situated  4.82  m 
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Figure  2.  Experimental  configuration  for  tomography  transmission  test  (Norsk,  1986). 


Table  2:  LOCATIONS  OF  ACOUCTIC  ELEMENTS 

Mooring 

Type 

Lat(®N) 

Lon(®E) 

A 

Source 

74.63 

23.40 

B 

Hydrophone  Array 

74.32 

23.55 

C 

Transceiver 

74.08 

22.00 

D 


Transceiver 


74.55 


21.72 


above  the  sea  floor.  Both  rav  and  mode  propagation  models  are  required  for  the 
identification  of  the  plane  wave  and  modal  beamfoimed  acoustic  arrivals  and  the 
subsequent  inversion  of  the  travel  time  data.  This  study  focuses  on  the  use  of  a  three- 
dimensional  ray  propagation  OKxlel  for  the  predictitm  of  the  beamformed  ray  arrival 
structure.  The  modeling  of  the  modal  arrival  structure  is  currently  in  progress  (Chiu,  et  al., 
1993(a)).  The  beamfmning  techniques  used  are  discussed  in  Miller,  et  al.,  (1993),  and  the 
hybrid  ray-noode  inversion  method  used  is  outlined  in  Chiu,  et  al.,  (1993(a)). 


C.  THESIS  OBJECTIVES  AND  APPROACHES 

The  objective  of  this  thesis  is  to  examine  and  establish  the  ray  forward  problem 
associated  with  the  BSPFEX  Coastal  TotxK>graphy  Test  A  three  dimensional  acoustic  ray 
theory  approach  was  used  to  simulate  the  transmission  of  the  224  Hz  tomographic  pulse 
signal  from  tfie  near  bottom  sound  source  (A)  to  the  vertical  receiver  array  (B), 
approximately  34  km  downrange.  To  dcvelc^  a  prediction  of  the  plane  wave  beamformed 
arrival  structure  a  dense  set  of  rays  was  first  traced.  Based  tm  these  rays,  eigenray  s  and  their 
anoplitudes  and  phases  were  then  determined.  The  predicted  arrival  structure,  as  a  function 
of  time-  and  angle,  is  the  key  for  identifying  the  observed  arrivals.  The  simulation  was  done 
fOT  both  a  three-dimensional  (3-D)  and  a  two-dimensional  (2-D)  modeled  ocean 
environment  in  order  to  examine  the  importance  of  the  3-D  effects. 

The  model  used  for  the  ray  tracing  was  the  Hamiltonian  Ray  Tracing  Program  for  the 
Ocean  (HARPO).  The  input  ocean  environments  were  mathematical  models  based  on  a 
bathymetric  chart  oi  the  region  (Norsk,  1986)  and  derived  sound  ^}eed  values  from 
seventy-two  evenly  spaced  Conductivity/Tempcraturc/Depth  (CTD)  stations  within  the  test 
area.  The  estimation  of  eigenrays  and  the  associated  arrival  structure  for  both  the  2-D  and 
3-D  environments  was  accomplished  using  a  MATLAB  program  that  operates  on  the  rays 
traced  by  HARPO. 
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D.  THESIS  OUTLINE 

The  remainder  of  this  thesis  consists  of  four  chapters.  Chapter  n  discusses  acoustic  ray 
theory,  the  ray  tracing  program  HARPO,  and  a  description  of  the  methods  used  to  calculate 
the  amplitude  and  phase  of  the  individual  arrivals. 

In  Chapter  m,  the  methods  used  to  tiKxlel  the  environments  far  input  into  HARPO  are 
described.  Chapter  IV  presents  and  discusses  the  simulated  arrival  stmcture  and  Chapter  V 
presents  the  conclusion  of  this  study. 
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a  ACOUSTIC  RAY  THEORY  MODELING 


A.  INTRODUCTION 

The  model  used  for  the  numerical  raytradng  was  the  NOAA  Hamiltonian  Acoustic 
Rayuracing  Program  for  the  Ocean  (HARPO)  (Jtmes,  et  al.,  1986).  HARPO  is  a  (3-D) 
FORTRAN  code  that  determines  acoustic  path  solutions  by  integrating  Hamilton's 
equation.  Jones,  et  al.  (1986)  have  provided  a  conq>lete  description  of  the  algorithm.  For 
input,  HARPO  requires  a  continuous  representation  of  a  3-D  sound  speed  field,  and  a 
continuous  2-D  representation  of  the  upper  and  lower  reflecting  surfaces,  Le.,  the  sea 
surface  and  bathymetry.  Gridded  bathymetric  and  sound  speed  data  can  be  made 
continuous  by  FORTRAN  subroutines  external  to  HARPO  through  the  use  of  empirical 
orthogonal  functions  (EOFs)  and  splines  (Qiiu,  et  al.,  1993(b)).  HARPO  itself  does  not 
provide  a  prediction  of  the  multipath  arrival  structure,  but  a  MATLAB  program 
'raystructure.m'  developed  by  Chiu  (1993)  can  be  applied  to  HARPO  output  files  for  the 
searches  of  eigenrays  and  the  computation  of  ray  amplitudes  and  phases  based  on  estimates 
of  ray  tube  spreading  and  losses  and  phase  modifications  due  to  boundary  interactions. 

B.  HAMILTONIAN  RAY  TRACING 

Acoustic  energy  travels  in  the  form  of  compressional  waves,  but  in  the  high  fiequency 
approximation  these  waves  behave  like  particles  whose  paths  are  described  by  integrating 
a  differential  form  of  the  wave  equatiotL  HARPO  bonows  Hamilten's  equation  from 
quantum  mechanics  and  q)plies  it  to  sound  propagation.  The  general  form  of  Hamilton's 
equation  is: 
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^  dH 
dt 

^  dH 
dt  5^ 


i=  1,2.3 
i  =  1.2.3 


(1) 


where  the  Hanultonian  function  H(pjfi2>P3:<Ij/l2>R3)  rq>resents  a  system's  total  energy 
with  respect  to  a  generalized  coendinate  system  p.  and  noomenta  q.  The  coordinate  system 
used  by  HARPO  is  spherical  polar,  and  fca*  acoustic  ^licatkHis.  q  represents  the  wave 
number  vector.  Therefore,  upon  specifying  initial  values  for  Equation  (1).  solutions  are 
obtained  by  integrating  the  system  of  differential  equadmis.  The  form  of  the  Hamiltonian 
in  sound  propagation  is: 

HiPi^qi)  =  (i>^-c(Pi)q^  (2) 

where  (O  represents  the  angular  wave  frequency.  c(Pi)  is  the  sound  speed  field,  and  ^  is  the 
magnitude  squared  of  the  wave  number  vector  (Jones  et  al..  1986). 


C.  ARRIVAL  STRUCTURE  ESTIMATION 

Within  the  ray  theory  framework,  the  received  signal,  rfr),  at  a  single  phone  can  be 
approximated  as: 

N 

r  (r)  =  OjS  (r  -  tj)  exp  (3) 

y  =  i 

where  N  is  the  number  of  eigenrays,  s(t)  is  the  signal  emitted  by  the  source,  and  a^-,  and 
tj  respectively  represent  the  anq)litude,  phase  shift  and  travel  time  for  each  eigenray  arrival. 

The  travel  time,  or  path  integral  of  sound  slowness,  was  calculated  in  HARPO  as  each 
ray  was  traced.  The  MATLAB  program*  raystracture.m’  was  used  to  find  eigenrays  and 
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determine  their  respective  an^litudes  and  phases  in  order  to  develop  an  estimate  of  the  ray 
arrival  structure  using  Equation  (3). 

F<nr  each  eigenray,  the  phase  shift  was  modified  for  each  boundary  interaction  and 
turning  point  (caustic  interaction).  For  each  surface  reflection,  the  phase  was  reduced  by  n . 
For  a  turning  point  (caustic  interacion),  it  was  reduced  by  7c/2,  and  for  a  bottom 
interaction,  the  phase  of  the  reflection  coefficient  evaluated  at  the  ray  angle  incident  on  the 
bottom  was  added. 

Eigenray  amplitude  was  calculated  by  ccnrecting  for  losses  due  to  ray  tube  spreading, 
surface  scattering  and  bottom  reflection.  A  description  of  these  factors  and  the  equations 
used  fcH-  the  correction  are  given  below. 


1.  Ray  Tube  Spreading 


Spreading  loss  is  the  geometrical  effect  that  represents  the  weakening  of  a  sound 
signal  as  it  spreads  outward  from  a  source  (Urick,  1983).  This  weakening  of  the  sigrud  is 
the  result  of  the  focusing  or  defocusing  of  adjacent  rays  as  they  propagate  down  range.  The 
transmission  equation  for  ray  tube  spreading  fcnr  a  non-absOTbing  medium  (Clay  and 
Medwin,  1977)  is: 


..2  _ 

■  o^OrAcIcosei 


(4) 


where  represents  a  reference  radius,  the  pressure  at  ,  Cq  is  tire  sound  speed  at 
and  6q  is  the  grazing  elevation  angle  at  Rf^  P,  c,  and  0  respectively  represent  the 
pressure,  sound  speed  and  angle  at  range  r.  The  vertical  separation  between  adjacent  rays 
at  range  r  is  represented  by  h,  and  A6  is  the  increment  between  the  irutial  launch  angles  of 
the  adjacent  rays.  The  MATLAB  program  ’raystructure.m'  uses  this  equation  to  calculate 
the  effea  of  ray  tube  spreading  on  relative  magnitude. 
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2.  Surface  Scattering 

The  shallowness  of  the  Barents  Sea  plays  a  significant  role  in  sound  propagation 
in  that  many  rays  interact  with  the  surface.  The  sea  surface  can  both  reflect  and  scatter 
sound.  Fot  a  smooth  sea,  the  magnitude  of  the  reflection  coefficient  (ratio  of  reflected 
pressure  to  incident  pressure)  is  unity.  For  a  rough  sea  surface  however,  because  some  of 
the  energy  is  scattered,  the  magnitude  of  the  coherent  reflection  coefficient  is  less  than 
unity  and  decreases  with  increasing  wave  height  (Brekhovsldkh,  et  al.,1982).  This  coherent 
reflectimi  coefficient  is  dependent  on  the  Rayleigh  parameter 

R  =  ikasinO  (5) 

where  k  is  the  wave  number,  a  represents  the  root  mean  square  displacement  of  a  rough 
surface,  and  6  is  the  grazing  angle  relative  to  the  horizontal.  The  magnitude  of  the  coherent 
reflection  coefficient,  for  an  irregular  surface  can  be  expressed  as  (Brekhovsldkh,  et 
al.,1982): 


R^fc-  expi-lR^)  (6) 

The  method  described  above  was  used  in  'raystructure.m'  to  include  surface  scattering 
losses. 

3.  Bottom  Loss 

Bottom  interaction  can  result  in  a  loss  of  acoustic  energy  due  to  both  scattering  of 
energy  out  of  the  ray  tube  and  the  transmission  of  some  of  the  energy  into  the  bottom.  This 
energy  loss  can  be  parameterized  by  the  reflection  coefficient  at  the  bottom.  If  a  plane  wave 
is  incident  on  a  boundary  of  two  fluids  with  density  and  and  sound  speed  Cj  and  C2, 

at  an  angle  of  6^ ,  then  the  snKwth  interface  reflection  coefficient  (Qay  and  Medwin,  1977) 
can  be  expressed  as: 

P2C2COS0j-PjCiCOS02  ^ 

~  P2C2COS0j  +  PjCiCOS02 
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where,  by  Snell’s  law: 


(8) 


The  program  'raystructiire.m'  calculates  the  bottom  reflection  coefficient  in  this  manner.  A 
bottom  sound  speed  of  3200  m/s  and  a  bottom  density  of  2700  kgAn^  were  used  in  this 
simulation.  Bottom  scattering  loss  was  assumed  to  be  negligible  and  therefore  not  modeled. 
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m.  MODELING  THE  BARENTS  SEA  ACOUSTIC  ENVIRONMENT 


A.  GENERAL 

Environmental  input  noodels  for  the  simulation  of  the  acoustic  ray  arrival  structure 
were  developed  using  data  that  provided  an  adequate  description  of  the  bathymetry  and 
sound  speed  fields.  The  accuracy  of  the  input  models  used  to  describe  the  ocean 
environment  determine  the  accuracy  in  tl:e  calculation  of  sound  propagation  solutions.  A 
physical  oceanognq)hic  descriptitm  of  the  experimental  region  is  given  in  Chapter  L  In  this 
chiq)ter,  the  construction  of  the  mathematical  environmental  models  input  into  HARPO  is 
discussed. 

B.  SOUND  SPEED  FIELD 

The  3-D  sound  speed  field  was  developed  using  measured  sound  speed  profiles  fiom 
72  CTD  casts.  The  locations  of  these  CID  stations,  occupied  during  the  first  ship  survey  in 
the  BSPFEX,  are  shown  in  Hgure  4.  Each  measured  profile  was  first  extended  to  a  uniform 
depth  and  then  interpolated  to  provide  data  at  a  two  meter  spacing  in  the  vertical  dirrction. 
Because  the  focus  of  this  study  is  on  the  propagation  path  between  the  224  Hz  source  and 

the  vertical  hydrophone  array,  the  modeled  region  was  chosen  to  be  a  36  x  20  krn^  area 
occupied  by  the  five  tracks  shown  in  Hgure  S.  Each  track  is  36  km  long  with  5  km 
horizontal  spacing  between  them.  Using  a  horizcmtal  inverse  distance  technique,  a  value  for 
sound  speed  was  interpolated  every  kilometer  along  each  track  for  every  two  meters  in 
depth.  The  resulting  profiles  were  further  srrx)othed  vertically  using  an  eleven  point 
hamming  filter.  This  is  required  to  eliminate  unrealistic  ray  chaos  resulting  from  the  break 
down  of  the  ray  approximation  in  the  presence  of  fine  stmcture.  Hgure  6  shows  the  sound 
q)eed  profiles  fiom  the  four  CTD  stations  located  along  the  line  of  sight  path  between  the 
source  and  the  array.  The  interpolated  and  smoothed  profiles  for  the  same  vertical  cross 
section  are  depicted  in  Figure  7.  A  contour  plot  of  the  sound  speed  field  along  the  track 
between  the  224  Hz  source  and  the  vertical  array  is  provided  in  Figure  8. 
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Figure  4.  Barents  Sea  Polar  Ront  Experiment  CTD  locations. 


Figure  5.  The  geographic  regicm  used  in  the  environmental  input  models  is 
rqnesented  by  the  five  tracks  shown  here.  Each  track  is  36  km  long  with  S  km 
separation  between  them. 
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Depth(m)  -  I  Deplh(fn) 


•'igure  6.  Sound  speed  profiles  firom  the  four  CTD  stations  along  the  line  of 
i^t  between  die  source  and  receiver. 


Figure  7.  Interpolated  and  filtered  sound  ^leed  profiles  along  the  line  of  sight 
b^een  the  source  and  receiver. 


IS 


Figure  8.  Sound  speed  ctmtour  alcmg  the  line  of  sight  between  the  source  and 
receiver.  (+1400  i^s) 


Figure  9  shows  the  sound  speed  contours  aloag  each  track.  Note,  from  Hgure  9,  that  to  the 
east  of  the  line  of  sight  track  lies  warmer  water  with  a  higher  sound  speed  which  will  tend 
to  re&act  sound  towards  die  west  Finally,  the  gridded  sound  speed  data  along  the  five 
tracks  was  made  continuous  for  input  into  HARPO  dirough  the  use  of  empirical  orthogonal 
fimctions  (EOFs)  and  splines  using  two  FORTRAN  subroutines  'gridder.f  and  'cgrid.f 
(Newhall,  et  al.,  1987).  To  determine  an  adequate  number  of  EOFs  needed  to  accurately 
represent  the  actual  sound  speed  field,  a  MAILAB  program  'eof m’  (Chiu,  1993)  was  used, 
calculating  the  maximum  error  associated  with  a  given  number  of  EOFs.  To  achieve  an 
acceptable  error  of  lO'^m/s,  sixteen  EOFs  woe  used  for  the  3'D  environment 
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Range  (km) 

Figure  9.  Sound  speed  contours  (+ 1400  m/s)  along  each  track  depicted  in 
Figure  5.  These  contours  represent  the  three  dimensional  sound  speed  field 
used  for  input  to  HARPO. 
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C.  BATHYMETRY  AND  SEA  SURFACE 

The  bathymeoy  input  fOT  HARPO  was  taken  from  the  chart  shown  in  Figure  2 
(Norsk,1986).  The  region  bounded  by  21.00  E,  24.00  E,  73.90  N  and  74.80  N  was  gridded 
into  76  subdivisions  on  each  side,  providing  a  0.45  km  resolution  in  the  longitudinal 
direction  and  0.9  km  resolution  in  the  latitudinal  directimt.  Water  depths  were  read  directly 
off  the  chart  at  those  grid  points.  Additionally,  Acoustic  Doppler  Current  Profiler  (ADCP) 
data  provided  depths  along  the  track.  These  two  sets  of  data  were  merged  to  provide  a  high 
resolution  bathymetry  model  for  the  region.  The  resulting  bathymetry  contours  are  shown 
in  Figure  10.  Because  the  modeled  envirmiment  provided  such  high  resolution,  it  was 


Figure  10.  Contour  of  bathymetry  model  used  as  input  for  HARPO  (depth  in  meters). 


necessary  to  smooth  it  using  a  low-pass  filter  to  eliminate  ray  chaos.  Similar  to  the  case  of 
the  sound  speed  field,  the  input  bathymetry  field  was  reduced  to  a  size  of  36  x  20  km^  and 
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regridded  to  one  kilometer  lesolutitMi  in  the  latitudinal  direction  and  five  kilometer 
resolution  in  the  Imigitudinal  direction  using  linear  interpolation.  A  3-D  view  of  the  input 
bathymetry  model  is  shown  in  Figure  11.  Superimposed  on  Hgure  1 1  are  the  locations  of 
the  224  Hz  source  and  the  hydrophone  array  as  well  as  a  bottom  interacting  ray  path.  Note 
that  the  depth  increases  to  the  southeast  which  will  tend  to  steer  bottom  interacting  rays  in 
that  direction,  counteracting  the  westward  refracting  effect  of  the  sound  speed  field.  Using 
two  other  FORTRAN  subroutines  'tgridder.f  and  ’bottom.f  (Newhall,  et  al.,  1987),  this 
gridded  input  was  splined  to  provide  a  continuous  bathymetry  representation  for  use  in 
HARPO. 


Figure  11.  Three  dimensional  representation  of  the  region  of  study.  Superimposed 
on  this  figure  are  the  locations  of  the  224  Hz  source  and  the  hydrc^hone  array  as 
well  as  a  bottom  interacting  ray  path. 
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regridded  to  one  kilcmieter  resolution  in  the  latitudinal  direction  and  five  kilometer 
resolution  in  the  longitudinal  direction  using  linear  inteipolation.  A  3-D  view  of  the  input 
bathymetry  model  is  shown  in  Figure  1 1.  Superimposed  on  Figure  1 1  are  the  locations  of 
tile  224  Hz  source  and  the  hydrophone  array  as  well  as  a  bottom  interacting  ray  path.  Note 
that  the  depth  increases  to  the  southeast  which  will  tend  to  steer  bottom  interacting  rays  in 
that  direction,  counteracting  the  westward  refracting  effect  of  the  sound  speed  field.  Using 
two  othw  FORTRAN  subroutines  'tgridder.f  and  'bottomi*  G^ewhall,  et  al.,  1987),  this 
gridded  input  was  splined  to  provide  a  continuous  bathymetry  representation  for  use  in 
HARPO. 


Figure  11.  Three  dimensional  representation  of  the  region  of  study.  Superimposed 
on  this  figure  are  the  locaticms  of  the  224  Hz  source  and  the  hydrophone  array  as 
well  as  a  bottom  interacting  ray  path. 


The  sea  surface  W9S  modeled  as  a  sidieie  concentric  widi  die  earth  with  a  radius  of 
6360.481  kan. 

D.  TWO  DIMENSIONAL  ENVIRONMENT 

The  2-D  environment  was  created  by  substituting  the  line  of  sight  bathymetry  and 
sound  speed  profiles  into  each  of  the  five  tracks.  Because  a  smaller  number  of  initial  sound 
speed  profiles  were  used  to  create  this  environment,  (mly  six  EOFs  were  necessary  to 
achieve  the  same  accuracy  as  the  3-D  case. 
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IV.  SIMULATION  RESULTS  AND  ANALYSIS 


A.  INTRODUCTION 

In  simnl^fing  the  propagation  of  die  224  Hz  tomognqihic  signal  from  the  near  bottom 
source  throu^  the  Barents  Sea  Ptdar  Fnmt  to  die  vertical  hydrophone  array,  rays  were 

traced  for  launch  angles  between  0.00^  and  40.00^  at  O.Q2S^  increments  for  both  the  2-D 
and  3-D  modeled  envircximents.  Eigenrays  and  their  amplitudes,  phase  shifts  and  travel 
times  were  then  estimated  for  ccmstructing  die  arrival  structure. 

B.  RAY  PATH  GEOMETRY 

All  rays  were  launched  along  a  constant  azimuth  and  traced  to  a  range  of 34.32493  km. 
Wavefront  corrections  were  applied  to  ray  path  length,  arrival  time  and  arrival  depth. 
Because  of  the  downward  refracting  sound  speed  profile  (Figure  6),  many  of  the  rays 
launched  at  shallow  and  intermediate  graang  angles  were  refracted-bottom-reflected 

(RBR).  Figure  12  shows  the  rays  launched  at  S.(X)^  for  bodi  the  2-D  and  3-D  cases.  In  both 

cases,  the  first  ray  to  interact  with  the  surface  was  launched  at  10.^5°.  These  rays  are 
dqiicted  in  Hgure  13.  The  first  ray  to  be  ctHupletely  suifoce-reflected  and  bottom-reflected 

was  the  13.35°  ray  for  both  the  2-D  and  3-D  cases.  Hgure  14  shows  these  rays.  The  3-D 
environment  introduced  horizontal  refraction  into  ray  path  geometries.  As  expected,  those 
rays  launched  at  shallow  angles  had  minimal  bottom  intnactions  and  were  driven  mainly 
by  the  sound  speed  field.  Due  to  the  warmer  water  to  the  east,  these  rays  were  refracted 
westward.  Those  rays  launched  at  steeper  grazing  angles  underwent  numerous  bottom 
intoactions  and  were  steered  eastward  by  the  bathymetry.  Figure  IS  shows  the  horizontal 
refraction  of  selected  rays. 
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Figure  13.  Path  geometiy  of  the  ray  launched  at  10.05°  for  (a)  the  2-D  environment 
and  (b)  the  3-D  environment 


(a) 


Figure  14.  Path  geometiy  of  the  ray  launched  at  13.35°  for  the  (a)  2-D  enviionmoit 
and  (b)  the  3-D  environment 
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but  clearly  shows  that  rays  launched  at  shallow  grazing  angles  are  driven  westwaid 
by  the  sound  speed  field  viieieas  those  rays  launched  at  steq>er  angles  are  steered 
eastward  by  the  bathymetry. 

C.  ARRIVAL  STRUCTURE 

To  examine  die  arrival  structure,  various  plots  were  made  based  on  all  of  the  rays 
traced.  Hgure  16  is  a  plot  of  arrival  angle  versus  arrival  depth.  Both  the  2-D  and  3-D  cases 
show  similar  forms  in  that  the  majority  of  the  rays  arrive  at  dq>ths  greater  than  100  m. 
Those  rays  that  are  not  completely  surface-reflected  and  bottom-reflected,  i.e.,  those  rays 

with  launch  angles  less  than  13.35°,  show  a  distinctive  relationship  between  launch  angle 
and  arrival  depth  in  that  the  minimum  depth  of  arrival  decreases  with  increasing  launch 
angle.  At  a  given  dq>th,  the  eigenrays  are  those  rays  launched  at  the  specific  angles 
determined  by  the  angle  versus  depth  curve.  In  general,  the  number  of  eigenrays  increases 
with  depth. 
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Figure  16.  Arrival  depth  versus  launch  angle  for  the  (a)  2-D  environment  and  the 
(b)  3-D  environment. 


26 


Figure  17  shows  a  plot  of  launch  angle  versus  arrival  time.  In  this  riguie  one  can  see 
that  the  earliest  arrivals  were  those  rays  launched  at  smaller  angles.  As  the  launch  angle 
increases,  so  does  the  cmresptxuiing  travel  time.  Note  that  fra*  the  3-D  case,  the  travel  tunes 
for  those  rays  launched  at  the  steepest  angles  are  a  little  less  than  those  launched  at  the  same 
angles  in  the  2-D  case.  This  is  likely  due  to  the  fact  that  the  rays  launched  at  steeper  angles 
undergo  more  boundary  interactions.  As  discussed  earlier,  the  bathymetry  in  the  3-D  case 
would  tend  to  steer  these  rays  towards  the  east,  a  region  where  sound  speed  is  higher.  From 

Figure  15,  the  ray  launched  at  40.00°  arrived  at  the  receiver  range  approximately  6°  east 
of  the  line  of  sight,  corresponding  to  a  distance  of  about  3.6  km  east  of  the  array. 

Figure  18  shows  a  plot  of  arrival  depth  versus  arrival  time.  Arrivals  less  than  25  s 
(corresponding  to  shallow  and  moderate  launch  angles)  show  much  more  regularity  than 
those  later  arrivals  associated  with  steeper  rays.  The  irregularities  among  the  arrivals 
greater  than  25  s  indicate  that  the  wave  front  has  become  distorted  as  a  result  of  the 
numerous  bottom  interactions.  The  implications  of  this  particular  effect  on  the  signal 
processing  beamfmming  will  be  discussed  later. 


C.  PLANE  WAVE  BEAMFORMED  ARRIVAL  STRUCTURE 

This  objective  of  tins  study  is  to  provide  a  prediction  of  the  plane  wave  beamformed 
arrival  structure  for  the  Barents  Sea  tomography  transmission  test.  This  prediction  can  be 
achieved  by  analyzing  the  arrival  angle  and  time  dependence  of  the  eigenrays  at  the  center 
of  the  array  (located  at  a  depth  of  approximately  195.8  m).  An  M.S.  thesis  study  by  Phil 
McLaughlin,  LT,  USN,  has  incorporated  these  predicted  arrivals  for  die  identification  of 
the  observed  plane  wave  beamfcHmed  arrivals. 

Figure  19  shows  the  arrival  structures  for  both  the  2-D  and  3-D  cases.  There  were  162 
eigenrays  identified  in  the  2-D  simulation  and  143  eigenrays  in  the  3-D  case.  A  comparison 
of  the  arrival  structures  for  the  two  simulations  is  given  in  Figure  20.  It  is  interesting  to  note 
that  for  arrivals  less  than  25  s,  both  cases  compare  vay  well  with  each  other.  Again,  these 
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Figure  18.  Arrival  depth  versus  arrival  time  for  (a)  propagation  in  a  2-D 
environment,  and  (b)  propagation  in  a  3-D  environment 
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Figure  19.  Channel  7  arrival  angle  versus  arrival  time  for  (a)  propagation  in  a 
2-D  environment,  and  (b)  propagation  in  a  3-D  environment. 
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arrivals  are  associated  with  rays  launched  at  low  and  moderate  angles.  After  25  s,  however, 
the  differences  between  die  two  cases  gradually  increase  and  the  latest  arrivals  differ  as 
much  as  300  ms.  These  latter  arrivals  are  associated  with  the  steep  ingle  rays  which, 
through  bottom  interaction,  undergo  a  significant  amount  of  horizontal  refraction.  The 
result  of  the  observed  plane  wave  plane  wave  beamformed  arrival  structure  is  given  in 
Figure  21  (Miller,  et  al.,  1993).  This  figure  shows  the  normalized  arrival  structure  for  plane 
wave  beamformed  data  at  one  degree  increments  from  -19.00^  to  19.00°.  The  dark  areas 
represent  the  peaks  of  the  arrivals.  Note  that  this  data  compares  favorably  with  die 
predicted  arrival  structure.  Figure  22  shows  a  plot  of  relative  magnitude  versus  arrival  time 
and  arrival  angle.  Both  cases  show  certain  similarities  in  that  the  relative  magnitude 
decreases  with  increasing  arrival  time.  This  dependence  is  primarily  the  result  of  surface 
scattering,  and  again,  because  the  later  arrivals  are  associated  with  steeper  angles,  they  are 
thus  subject  to  more  boundary  interactions,  hence,  greater  losses.  Based  on  the  work  by 
Emblidge  (1991),  mean  wave  heights  of  one  meter  were  used  in  modeling  sea  surface 
roughness. 
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Figure  20.  Channel  7  predicted  plane  wave  beamformed  arrival  structure.  (2-D  results 
are  represented  by  "o”,  3-D  results  are  represented  by  "+"). 
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Figure  21.  Coherent  average  of  28  sequences  of  len^  3.9375  seconds  for  plane 
wave  beainforming  angles  from  -19  to  19  degrees  in  increments  of  one  degree  (from 
Miller,  et  al.,  1993). 
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Figure  22.  Relative  magnitude  of  Qiannel  7  airivals  for  (a)  propagation  in  a  2-D 
environment,  (b)  propagation  in  a  3-D  environment 


34 


Arrival  Angle  (deg) 


V.  CONCLUSIONS 


It  was  the  objective  of  this  thesis  to  en^loy  a  ray  theory  model  to  produce  the 
predicted  plane  wave  beamfoimed  arrival  structure  £»-  the  1992  Barents  Sea  tomography 
transmissitHi  test  Bodi  a  2-D  and  3-D  modeled  environment  were  used  in  an  effort  to 
determine  the  effects  of  horiztmtal  refractitm  (»i  the  arrival  structure. 

For  both  cases,  arrivals  could  be  separated  into  two  distinct  groups: 

1)  Early  Arrivals  -  those  arrivals  with  travel  times  less  than  or  equal  to  25  s.  These  arrivals 
are  associated  with  rays  launched  between  0.00^  and  28.00°  and  have  arrival  angles 
between  -14.00°  and  14.00^. 

2)  Late  Arrivals  -  those  arrivals  with  travel  times  greater  dran  25  s.  These  arrivals  are 
associated  with  launch  angles  greater  than  28.00^  and  have  arrival  angles  between 
±15.00‘’and  ±20.00^*. 

In  both  the  2-D  and  3-D  cases,  the  late  arrivals  showed  a  distortion  of  the  wave  front 
In  terms  of  signal  processing  inq)lications,  this  could  reduce  the  accuracy  with  which  the 
travel  time  is  estimated  for  the  steep  angle  rays  using  plane  wave  beariiformmg. 

As  ejqrected,  in  both  cases,  later  arrivals  exhibit  significantly  more  loss  than  early 
arrivals  due  to  numerous  boundary  interactions. 

The  effea  of  prc^gation  in  the  3-D  environment  appears  to  be  prominent  only  in  the 
late  arrivals.  This  is  primarily  due  to  the  nature  of  the  bottom  which  is  causing  both  wave 
front  distortion  and  significant  horizontal  refraction  for  the  very  steep  angle  rays. 

Most  inq)ortantly  however,  this  study  has  denx>nstrated  that  an  acoustic  ray  theory 
model  can  be  used  to  accurately  predict  the  plane  wave  beamformed  arrival  structure  for 
shallow  water. 


35 


LIST  OF  REFERENCES 


Barents  Sea  Polar  Front  Group,  Experimental  Plan  Barents  Sea  Polar  Front  Experiment 
August  1992,  May  195^.  (Bouike,  RJL,  C-S.  Chiu,  J.F.  Lynch,  JJi.  Miller,  RD. 
Muench). 

Brekhovskikh,  L.  and  Yu.  Lysanov,  Fundamentals  of  Ocean  Acoustics.  Springer- Verlag 
Berlin  Hiedelbuig,  New  York, 1982. 

Clay,  CS.  and  H.  Medwin,  Acoustical  Oceanography,  Wiley  &  Sons,  New  York, 
19T7. 

Chiu,  C.-S.,  Department  of  Oceanography,  Naval  Postgraduate  School,  Monterey,  CA, 
Personal  Communicadon,  1992-1993 

Chiu,  C.-S.,  JJL  Miller,  and  J.F.  Lynch,  “Inverse  Techniques  for  Coastal  Acoustic 
Tomography”,  submitted  to  the  Proceedings  of  the  1993  International  Conference 
on  Theoretical  and  Con[q)utational  Acoustics,  July  1993. 

Chiu,  C.-S.,  A.J.  Semmer,  C.M.  Ort,  JJL  Miller  and  LL.  Ehret,  “A  ray  variability  analysis 
of  sound  transmission  from  Heard  Island  to  California”,  accepted  by  the  Journal  of 
the  Acoustical  Society  of  America:  The  Heard  Island  Experiment  Special  Issue, 
1993. 

Elliott,  J.M.,  Simulation  of  Acoustic  Multipath  Arrival  Structure  in  the  Barents  Sea,  Naval 
Postgraduate  School  Thesis,  1992. 

Emblidge,  J.M.,  A  Feasibility  Study  of  Ocean  Acoustic  Tomography  in  the  Barents  Sea, 
Naval  Postgraduate  School  Thesis,  1991. 

Jones,  R.M.,  J.P.  Riley  and  T.M.  Georges,  HARPO  a  versatile  three-dimensional 

hamiltionian  ray-tracing  program  for  acoustic  waves  in  an  ocean  with  an  irregular 
bottom.  Wave  Propagation  Laboratoty,  NOAA,  Boulder,  Colorado,  457  pp.,  1986. 

Leong,  R,  Features  of  the  physical  oceanographic  conditions  of  the  Barents  Sea, 
unpublished  manuscript,  1991. 

Miller,  JJL,  Department  of  Electrical  and  Computer  Engineering,  Naval  Postgraduate 
School,  Mcmterey,  CA,  personal  communication,  1992-1993. 


36 


Miller,  Jii,  C-S.  Chiu  and  J.  F.  Lynch,  “Signal  Processing  far  Coastal  Acoustic 
T(»Dography”,  submitted  to  the  Proceedings  of  the  1993  Ibtemational  Conference 
cm  Theoretical  and  Computational  Acoustics,  August  1993. 

Munk  W.,  and  C  Wunsch,  Octcat  acoustic  tomography:  a  scheme  for  large  scale 
monitorings  Deep  Sea  Res.,  v.  26,  Part  A,  123-161, 1979. 

Newhall,  A£.,  JJ.  Lynch,  C-S.  Chiu  and  J.R.  Daugherty,  “In^novemcnts  in  the  threc- 
dimensitHial  raytracing  codes  for  underwater  acoustics.  Computational  Acoustics, 
V.  1, 1987. 

Norsk  Polarinstitutt,  Bathymetry  Barents  Sea,  Chan  No.  7421,  Olso,  Norway,  1986. 

Urick,  Robert  J.,  Principles  of  Underwater  Sound,  3rd  ed.  McGraw-Hill  Book  Co.,  New 
York,  1983. 


37 


INITIAL  OISTRlBUnON  LIST 


Defense  Technical  Infiormatira  Center 
Cameron  Station 
Alexandria,  VA  22304-6145 

Librarian 
Code  52 

Naval  Postgraduate  School 
411  Dyer  Road 
Monterey,  CA  93943-5002 

Oceaiiognq;>hy  Department 
CodeOaCo 

Naval  Postgraduate  School 
833  Dyer  Road 
Monterey,  CA  93943-5114 

Professor  Ching-Sang  Chiu 
CodeOOa 

Naval  Postgraduate  School 
833  Dyer  Road 
Monterey,  CA  93943-5 1 14 

Professor  James  R  Miller 
Code  EC/Mr 

Naval  Postgraduate  School 
833  Dyer  Road 
Monterey,  CA  93943-51 14 

Office  of  the  Oceanogrs^her  of  the  Navy 

Attn:  Dr.  Robert  ^^nokur 

Naval  Observatory 

34th  and  Massachusetts  Avenue  NW 

Washingtcm,  DC  20390 

Commander 

Naval  Oceanogrqrhic  Command 

Stennis  Space  Center 

Bay  St  Louis,  MS  39539-5000 


1 


8.  Dr.  Thomas  Cuttin 
Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  VA  22217-5000 

9.  LT  John  L.  Myl^  1 

Fle^  Numerical  Oceanography  Center 

7  Grace  Hopper  Avenue,  Stop  1 
Mcmterey,  CA  93943-5501 

10.  Professor  James  F.  Lynch  1 

Departmoit  of  Applied  Ocean  I%ysics  and  Engineering 

Woods  Hole  Oceanographic  Institution 
Woods  Hole,  MA  02543 

11.  Nfr.  Barry  P.  Blumentiial  1 

Office  of  Naval  Research 

800  North  Quincy  Street 
Arlington,  VA  22217-5660 

12.  Dr.  M(4isen  Badiey  1 

Office  of  Naval  Research 

800  North  Quincy  Street 
Arlington,  VA  22217-5000 


39 


